We study the H2O chemistry in star-forming environments 
under the influence of a central X-ray source and a central far 
ultraviolet (FUV) radiation field. The X-ray models are ap- 
phed to envelopes around low-mass Class and I young stel- 
lar objects (YSOs). The gas-phase water chemistry is mod- 
eled as a function of time, hydrogen density and X-ray flux. 
To cover a wide range of physical environments, densities be- 
tween riu = lO'^-lO^cm"^ and temperatures between T = 10- 
lOOOK are studied. Three different regimes are found: For 
T < 100 K, the water abundance is of order 10^^-10"^ and 
can be somewhat enhanced or reduced due to X-rays, depend- 
ing on time and density. For lOOK < T < 250 K, H2O is 
reduced from initial x(H20) ~ 10 following ice evapora- 
tion to x(H20) ~ 10"'' for Fx > 10"^ ergs s"' cm"^ (f = 
10"* yrs) and for Fx ^ 10""* ergs s"' cm^^ ^ lO^yrs). At 
^) , higher temperatures (T > 250 K) and hydrogen densities, wa- 
PsJ ' ter can persist with x(H20) ~ 10 even for high X-ray fluxes. 

Water is destroyed in both Class and I envelopes on rela- 
(]J tively short timescales (f » 5000 yrs) for realistic X-ray fluxes, 
, although the effect is less prominent in Class envelopes 
\^ due to the higher X-ray absorbing densities there. FUV pho- 
, tons from the central source are not effective in destroying 
water. X-rays reduce the water abundances especially in re- 
gions where the gas temperature is T < 250-300 K for fluxes 
Fx S lO^^-lO""* ergs s"' cm"^. The affected regions can be en- 
, velopes, disks or outflow hot spots. The average water abun- 
dance in Class I sources for Lx > 10^^ ergs s ' is predicted to 
O , be a:(H20) < 10 ^ Central UV fields have a neghgible influ- 
' ence, unless the photons can escape through cavities. 
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Abstract. 



1. Introduction 

The importance of water in star-formation is unquestioned. 
Water is found abundantly on dust grains as well as in the gas- 
phase of molecular clouds, cores, envelopes and protostellar 
disks. It is a key molecule in the hot (T > 100 K) gas-phase 
chemistry and is one of the main gas coolants. The evolution 
over time of molecular cores influences the abundance of H2O 
(Ceccarelli et al. 1996 1, making it sensitive to the chemical age 
of a source. Through its masing ability, water has been used to 
study dynamical properties of star-forming regions. 

In cold gas (T < 100 K), water is mainly produced by 
ion-molecule reactions. At temperatures above 250 K, H2O is 
largely produced in reactions of O and OH with H2. Water 
is also efficiently generated on dust grains where it forms an 
icy mantle. Observed abundances of water ice are typically 
x(U20) ^ 10 with respect to H2 (Tielens et al. 1991 ; Boogert 
et al. I2004t . At temperatures T > 100 K, H2O evaporates 
into the gas-phase and has a major impact on the chemistry 
(van Dishoeck & Blake Il998l, Doty et al. 2002 J . Water is also 
believed to be formed in shocked regions where it may be 
produced due to the high gas temperature and/or may be re- 
leased from grains. Towards high-mass star-forming regions, 
hot (T > 100 K) gas-phase water is usually observed with 
abundances jic(H20) = 10 comparable to the ice abundance 
- whereas at temperatures below lOOK, it is found to be ap- 
proximately 100 times less abundant (e.g., Boonman & van 
Dishoeck l2003l van der Tak et al. I2005t . The observed water 
emission towards low-mass objects has been interpreted with 
similar jumps at T a: 100 K (CeccarelH et al. 120001 Maret 
et al. I2OO2I 1. Although the solid water abundances are com- 
parable in low and high-mass regions, the hot gas-phase wa- 
ter in low-mass sources seems to be much less abundant than 
in the massive objects. In addition, a fraction of the emission 
is believed to arise in outflow regions, as inferred from ob- 
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servations both on and off source with the Long Wavelength 
Spectrometer (LWS) on board the Infrared Space Observatory 
(ISO) (Giannini et al. 2001. Benedettini et al. 2002 1. Nisini et 
al. J2002t compared the far-infrared spectra of both Class and 
Class I sources. Water lines were found to be prominent in the 
spectra of Class objects but they were not detected toward 
Class I sources corresponding to an upper limit on the abun- 
dances ofx(H20) < 10 ^ 

Protostars are often found to be sources of strong X-ray 
fluxes (e.g., Casanova et al. 1995 , Koyama et al. 1996 Imanishi 
et al. 2001 1. The observed X-rays are understood to be ther- 
mal emission from magnetic stellar activities or from the disk- 
star system (e.g., Feigelson & Montmerle 1999). Typical X- 
ray luminosities range from Lx ~ 10^^-10^' ergs s"' in the 
0.5-lOkeV band with plasma temperatures between 0.6-7 keV 
(e.g., Imanishi et al. 2001 1. Flares can even lead to plasma 
temperatures temporarily exceeding 10** K with luminosities 
higher than Lx ~ 10^^ ergs s '. Within the Chandra Orion 
Ultradeep Project (COUP), Preibisch et al. (2005 1 found a ra- 
tio for the X-ray luminosity to the bolometric luminosity of 
log(Lx/i'boi) ~ -3.6 for Class I and older protostars. No X- 
ray detection has been reported to date toward Class ob- 
jects (Hamaguchi et al. 2005a 1. However, this might be due to 
the large X-ray absorbing hydrogen and dust column densities 
found in these sources. The nature of the sources in recent re- 
ports of X-ray detections towards very young objects is poorly 
known and they could be either Class or I objects (Hamaguchi 
et al. 2005b Forbrich et aL l2005l . 

The influence of a central X-ray source on the envelope 
around young stellar objects has recently been studied by 
Stauber et al. ( 20fl5j . It was found that the ionization rate in 
the inner region of the surrounding envelope may be domi- 
nated by the ionizing X-ray flux rather than by an inner far 
ultraviolet (FUV) radiation field or the cosmic -ray ionization 
rate. The ionizing flux enhances H3 and He"^ that trigger a dis- 
tinct chemistry. It was also found that the H2O abundances 
in the gas-phase may be reduced due to X-rays. The abun- 
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dance and behavior of H2O in both the solid and the gas-phase 
will be studied extensively with the upcoming Herschel Space 
Observatory. It is therefore timely to investigate the influence 
of X-rays and inner FUV fields on the abundances of gas-phase 
water in more detail. 

In this paper, we first study the general influence of X- 
rays on gas-phase water. The time-dependent chemical X-ray 
models of Stauber et al. J2005> are used and applied to differ- 
ent density and temperature regimes for different X-ray fluxes 
(Sect. IS. The densities and temperatures are chosen to cover 
a wide range of conditions that are applicable to cold molecu- 
lar clouds, envelopes, protostellar disks and outflow hot spots. 
Regions with density and temperature gradients are modeled 
and discussed in Sect. |3 These regions are representative of 
low-mass Class and I envelopes. In addition to X-rays, cen- 
tral FUV fields are considered. The aim is to find systematic 
trends between the different type of sources. The results for the 
different envelopes are discussed in Sect. |3 The conclusions 
are drawn in Sect.|5l 

2. General parameter study 

The gas-phase water abundance is modeled as a function of 
time, hydrogen density (nn - n(H) + 2n(H2)) and X-ray flux. 
The results of this study can be applied to a wide range of dif- 
ferent physical environments including protostellar envelopes 
and disks. 

2.1. Chemical model 

The X-ray chemistry models are described in detail by Stauber 
et al. J2005t . The models are an extension of the gas-phase 
chemical models of Doty et al. ('2DD2"20D4^ to allow the impact 
of X-rays from the central source on the suiTounding molecu- 
lar envelope. The input parameters for the model are the initial 
molecular or atomic abundances, the hydrogen density, the gas 
temperature, the X-ray luminosity, the X-ray emitting plasma 
temperature, the X-ray absorbing hydrogen column density, the 
cosmic -ray ionization rate and the enhancement of the inner 
and outer FUV field with respect to the average interstellar ra- 
diation field (ISRF). 

The chemical model is based on the UMIST gas-phase 
chemical reaction network (Millar et al. ,1997j and calculates 
the time-dependent number density n{i) of each species for a 
certain temperature and distance from the source. For the ini- 
tial chemical abundances we follow the models of Doty et al. 
d2002ll2004l for flie high-mass source AFGL 2591 and low- 
mass source IRAS 16293-2422. These models successfully re- 
produced many of the observed molecular lines. The effects of 
evaporation of a certain species into the gas-phase have been 
approximated by initially depleting this species below its evap- 
oration temperature T^y. Specifically, all H2O is assumed to 
evaporate at T > lOOK. No photode sorption of ice is taken 
into account. The initial abundances are listed in Table The 
adopted cosmic-ray ionization rate is discussed in Sect. 13. II 

The gas temperature has been taken to be closely coupled 
to the dust temperature. Although the X-ray flux and the gas 
temperature are not fully independent, they have been treated 



as uncoupled variables in order to study the dependence of H2O 
on each of them. To cover a wide range of different physical en- 
vironments, the hydrogen density is varied between nu = 10"*- 
10^ cm"^. At higher densities, the mean free path of X-rays 
becomes small (a few AU). Typical densities in the inner part 
of envelopes or protoplanetary disk atmospheres are between 
hh = 10^-10^ cm"-', whereas those for outflows are rtu = 10"*- 
10^ cm"-'. The adopted gas temperature is between 10-1000 K. 
At higher temperatures (T > 2000 K), molecular hydrogen and 
water are dissociated due to H and H2 collisions. The modeled 
X-ray fluxes are between Fx = 10"''-10ergs s"' cm"^. The ini- 
tial X-ray spectrum is reduced preferentially at low energies 
by the intervening hydrogen column density such that the local 
X-ray flux roughly scales with 



\ lO^^i ergs s-'/\5.0x 1021 cm-2/ \56AU/ 

in ergs s"' cm"^, where Lx is the original X-ray luminosity, A^h 
the X-ray absorbing column density and r is the distance from 
the source. The X-ray fluxes are accurate within 20% for Tx - 
3 X 10^ K and A^H = 5 X lO^'-lO^^ cm-^ with y = 0.3, for A^h = 
2-5 X 10^3 cm-2 with y = 0.4 and for A^h = 6 x lO^^-lO^^ ^,^^-2 
with y - 0.5. A column density of A^h = 5 x 10^' cm"^ absorbs 
all photons below * 1 keV (Stauber et al. (20051 . The X-ray 
spectrum is assumed to be thermal with a plasma temperature 
Tx = 3 X 10^ K (2.6 keV). The shape of the X-ray spectrum, 
however, has only little influence on the chemistry (Maloney et 
al. rT996l Stauber et aL GOOSI l. FUV fields are neglected. 

All reaction rates that are relevant for the H2O chemistry 
(Sect. l23l are given with a maximum uncertainty of 25% by 
the UMIST database'. Self-shielding of CO has been included 
by using the shielding functions of Lee et al. (,1996il . 

2.2. Results 

The model results for different X-ray fluxes and hydrogen den- 
sities are shown in Figs. [2 and |2 for t = lO^yrs and in Figs.|5] 
and 0] for t - lO'^yrs. It can be seen in Figs. [0 and |5] that 
X-ray fluxes Fx ^ 10"^ ergs s"' cm"^ have no or minor in- 
fluence on the water abundances for the densities considered. 
For higher X-ray fluxes, three characteristic regimes can be dis- 
tinguished: 1. the regime at T < lOOK, where gas-phase wa- 
ter is mainly formed and destroyed in ion-molecule reactions 
(.^:(H20) ~ 10-^-10-^), 2. the regime 100 K < T < 250 K, 
where water is released from grains, but destroyed by X-rays 
(;c(H20) ~ 10"*) and 3. the regime T > 250 K, where the water 
abundance is x(H20) > 10 The three regimes are discussed 
in some detail in the following paragraphs. 

At t = lO'^yrs, X-ray fluxes 10"^ < Fx < 
10^2 ergs s"' cm"2 enhance the water abundance for tempera- 
tures T < 100 K compared to the model without X-rays. The 
reason for this is the increased HaO^ abundance due to X-rays. 
HaO^ can recombine quickly to H2O in reactions with elec- 
trons. Note that net formation of H2O at these fluxes can occur 
because H^O^ does not require H2O in its formation. It is pro- 
duced by ion-molecule reactions starting from CO, CO2, O, 

' http;//www.rate99.co.uk 
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Table 1. Initial abundances and cosmic-ray ionization rate. 



Species 


Initial abundance 


Ref. 


Initial abundances T > 100 K 






CO 


2.0E-04 


a 


CO2 


3.0E-05 


b 


H2O 


1.5E-04 


c 


H2S 


l.OE-08 


d 


H2CO 


8.0E-08 


d 


N2 


7.0E-05 


e 


CH4 


l.OE-07 


e 


C9H4 


8.0E-08 


e 


C2H6 


l.OE-08 


e 


CH3OH 


1.5E-07 


d 





0.0 


e 


s 


0.0 


e 


Initial abundances T < 100 K 






CO 


2.0E-04 


d 


CO2 


0.0 


f 


H2O 


0.0 


f 


H2S 


0.0 


f 


N2 


7.0E-05 


e 


CH4 


l.OE-07 


e 




8.0E-08 


e 


C2H6 


l.OE-08 


e 


H2CO (60 < r(K) < 100) 


8.0E-08 


d 


H2CO (r(K) < 60) 


0.0 


d 


CH3OH (60 < r(K) < 100) 


1.5E-07 


d 


CH3OH (r(K) < 60) 


0.0 


d 





l.OE-04 


d 


s 


6.0E-09 


g 


Cosmic-ray ionization rate 4i (s"') 


0.8E-17 


h 



All abundances are relative to molecular hydrogen. " l0rgensen et 
al. 1 2002 1, * Boonman et al. <2003at . Boonman & van Dishoeck 
(2003 I, Doty et al. l 2004"i, ' Chamley ('19971, f assumed to be 
frozen-out or absent in cold gas-phase, ' Doty et al. (.2002^ . '' see 
text. 



O2 and OH, which exist in the gas phase. Higher X-ray fluxes, 
however, can enhance water only at densities hh <: 10^ cm"^ 
for T < 100 K (Fig.|2j- At f = 10^ yrs, water is only enhanced 
at hh = 10^ cm"'' in this temperature regime (Fig.|5}. At higher 
densities, water is less abundant in models with X-rays com- 
pared to those without X-rays. 

For lOOK < T < 250 K, the initial abundance of water is 
much higher due to ice evaporation but this high abundance 
rapidly decreases for Fx ^ lO^-'ergss"' cm"^ at f = lO'^yrs 
and for Fx <: 10""^ ergs s"' cm"^ at f = 10^ yrs. H2O is mainly 
destroyed in reactions with H3 and HCO^ whose abundances 
are enhanced due to the X-rays. Another destruction mecha- 
nism is the dissociation of H2O by FUV photons created locally 
by excited H2. The rate for the X-ray induced FUV destruction 
of water is calculated with the method and numbers provided 
by Maloney et al. J1996> . It is an adaption of the treatment of 
Gredel et al. (1989'^ for the efi'ects of internally generated FUV 
photons by cosmic rays. 

In the temperature regime T > 250 K, the reaction OH 
H- H2 — > H2O + H can become more eflicient than the wa- 
ter destroying reactions and H2O has abundances of the order 



of lO"**. At low densities (hh = lO'^-lO^ cm"-'), high X-ray 
fluxes lead to the destruction of water even for high tempera- 
tures. Although the gas temperature at these high fluxes may be 
more than lOOOK, water does not reach fractional abundances 
higher than x(H20) » 10"^, since H2O is destroyed in colli- 
sions with H and H2 at higher temperatures. The water abun- 
dance is therefore x(H20) < 10"'" for Fx = 10 ergs s ' cm ^ 
and «H = lO'^cm"^ and ;c(H20) < 10"^ for «h = 105cm"^ 
independent of the gas temperature. At higher densities, how- 
ever, even high X-ray fluxes cannot destroy water in the regime 
T > 250 K. At densities «h S lO^cm"^, X-rays even enhance 
the water abundance (.ji:(H20) > 3 x lO""') for temperatures 
T > 300-600 K compared to models without X-rays. This en- 
hancement is again mainly through recombination reactions of 
H3O+. 

At later times (f = 10^ yrs), water is destroyed in the regime 
lOOK < T < 250K for Fx > 10"^ ergs s"' cm ^. The other 
temperature regimes in the X-ray models resemble those at f = 
10^ yrs. The main differences with the models at f = lO'^yrs, 
however, are the water abundances of the models without X- 
rays. H2O is gradually destroyed in reactions with HCO^ for 
temperatures T > 100 K. In models without X-rays, HCO^ is 
mainly formed by cosmic -ray induced reactions. At f = 10'' yrs, 
water between T x lOOK and T « 200 K is destroyed down 
to abundances of a few xlO"^. The water abundances for T < 
lOOK on the other hand, are slightly higher compared to the 
models at f = 10^ yrs. This is due to the relatively slow ion- 
molecule reactions that need time to build up water. 

2.3. Chemical reactions relevant for H2O 

Water is assumed to be initially frozen out on the dust grains for 
temperatures T < 100 K and injected into the gas at T = 100 K 
with a fractional abundance ji:(H20) = 1.5 x 10""*. In the pres- 
ence of X-rays, water is then mainly destroyed in reactions with 
H3 and HCO^. Another destruction route is the dissociation of 
water due to X-ray induced FUV photons. To estimate the im- 
portance of this reaction on the H2O abundance, models have 
been run for different X-ray luminosities ignoring this reaction. 
The resulting H2O abundances vary by less than 15% in com- 
parison to models where this reaction is included. It empha- 
sizes that the main destruction reactions of H2O are the ion- 
molecule reactions with HJ and HCO^, respectively. The bulk 
of the H2O not evaporated from the grains stems from the re- 
action of OH with H2, which becomes dominant for tempera- 
tures T > 250 K, and from electron recombination reactions of 
HaO^. The adopted branching ratios for the dissociative recom- 
bination of HsO^ to H2O and OH are 0.3 and 0.7, respectively 
(e.g., Jensen et al. 2000 1. H^O^ is abundantly produced in the 
reactions of water with H^ and HCO^ in the inner region of 
the envelope. Further away from the source, H30^ is most ef- 
ficiently formed in the following way: H^ reacts with atomic 
oxygen to produce OH^. This is followed by a series of hy- 
drogen abstract reactions with H2 leading to HaO^. The main 
source for H^ is the reaction of molecular hydrogen with Hj, 
which is in predominantly produced by electron impact ioniza- 
tion of H2. HCO^ is produced via the reaction of HJ with CO. 
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Fig. 1. Fractional water abundances as functions of the gas tem- 
perature and low X-ray flux (Fx - 10"''-10"'' ergs s"' cm"^) 
for different densities (cm^^) at f = lO^yrs. The vertical line 
indicates the 250 K temperature mark. 
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Fig. 2. Fractional water abundances as functions of the gas tem- 
perature and high X-ray flux (Fx = 10^^-10 ergs s"' cm"^) for 
diff'erent densities (cm^^) at f = 104yrs. The vertical line indi- 
cates the 250 K temperature mark. 



OH is largely produced in recombination reactions of H30"^ 
and destroyed in reactions with atomic oxygen or molecular 
hydrogen. The chemical network with reactions and species 
relevant for H2O in the X-ray models is presented in Fig. |5] 
Photoionizations and dissociations due to inner or outer FUV 
fields are not shown. 

3. Applications to protostellar envelopes 

In order to study the influence of X-rays on the gas-phase 
water abundance in more realistic environments, the models 
are applied to envelopes of Class I (t a 10^-10^ yrs) and 
Class (f ~ a few xl03-104yrs) low -mass YSOs. These re- 
gions may be characterized by a power-law density distribution 
and an increasing gas temperature towards the central proto- 
star. Class and I envelopes mainly diff'er in their total mass. 
Typical envelope masses of Class objects are Menv S 0.5 Mq 
whereas Class I sources have less massive envelopes with 
Menv ^ 0.5 Mq (e.g., J0rgensen et al. 'SUOT, '2005a'). Although 
spherical ID models are an approximation to the real geometry 
of protostellar cores, they successfully explain many of the ob- 
served features of low-mass protostars from few hundred AU 
out to 104 AU scales (e.g., Schoier et al. 2004,; J0rgensen et 
al. 120021 [2005c> and provide a framework for discussing their 
overall chemical structure (e.g.. Doty et al.'2004'; Maret et al. 
[2004.. J0rgensen et al. .2004j . Inside a few hundred AU, how- 



ever, the geometry is likely to be much more complex with cav- 
ities and a flattened protostellar disk in addition to the spheri- 
cal inner envelope. Since no realistic physical models based 
on observations on those scales are yet available, we adopt the 
spherical models extrapolated to small radii as a starting point. 
Possible efifects due to cavities and non-spherical symmetry are 
discussed throughout this section and in Sect.|3] 

In the following models, the H2O chemistry is assumed 
to be irradiated by X-rays and FUV photons from the cen- 
tral source. The FUV models have been described in detail by 
Stauber et al. t2004/ 

3.1. Class I envelope model 

The density and temperature profiles adopted for a prototypi- 
cal Class I source are presented in Fig.|6land Table |2l based 
on the observations by J0rgensen et al. (2 002) . In that anal- 
ysis, it was assumed that the gas temperature equals the dust 
temperature. X-rays and FUV fields, however, are able to in- 
crease the gas temperature due to photoelectric heating. The 
gas temperature has therefore been explicitly calculated for the 
different X-ray and FUV fluxes. The X-ray luminosity is varied 
between Lx = 10^^-10^" ergs s"' and the FUV field from the 
central source (Go in) at the inner radius from 10^ to 10^ with 
respect to the average ISRF. The heating rate due to X-rays 
is taken from Maloney et al. (iil996< . and for the FUV heat- 
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Fig. 3. Fractional water abundances as functions of the gas tem- 
perature and low X-ray flux (Fx - 10"''-10"'' ergs s"' cm"^) 
for different densities (cm^^) at f = 10^ yrs. The vertical line 
indicates the 250 K temperature mark. 
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Fig. 4. Fractional water abundances as functions of the gas tem- 
perature and high X-ray flux (Fx = 10^^-10 ergs s"' cm"^) for 
diff'erent densities (cm^^) at f = 10^ yrs. The vertical line indi- 
cates the 250 K temperature mark. 



ing we follow Bakes & Tielens f 1994V The heating rates are 
then compared to the cooling rate due to gas-dust collisions 
provided by Hollenbach & McKee (1989 1. At lower densities 
(riii X 10^ cm"^) line cooling due to C"^, O and CO also be- 
comes important (see also Doty & Neufeld ,1997 1, but these 
effects are minor for our purpose since the densities in the re- 
gions of interest are always «h > 10^ cm"-^. 

It is found that only X-ray luminosities Lx > 10^' ergs s"' 
lead to a noticeable increase in the gas temperature (Tgas - 
Tdust ^ 10 K). The results for Lx - 10^' ergs s ' and Lx = 
lO^^ergss^i are plotted in Fig.|6l The temperature of the gas 
is only affected in the inner » 0.1 AU by the central FUV field 
(not shown in Fig.|6|l. High FUV fields (Go.in <: 10^) are needed 
though to increase the gas temperature significantly above the 
dust temperature. 

The effects of the outer FUV field can be neglected. The 
optical depth in the region of interest (T > 100 K) is Ay ~ 18 
and all FUV photons should be absorbed by the outer envelope 
for typical values of the ISRF. The outer FUV field is taken to 
be G(),out = 1 according to the average ISRF. 

The adopted value of the cosmic-ray ionization rate ^cr = 
0.8 X 10"'^ s^i is based on the results of Stauber et al. ("2005') 
who showed that low- J lines of HCO^ tend to trace the cosmic- 
ray ionization rate, whereas high-7 lines are more sensitive to 
the ionizing inner X-ray flux. To constrain ^ci , the rate is var- 
ied in our models in order to fit the HCO"^ 1-0 observations 



He+ 




Fig. 5. Most important reactions involved in the X-ray chem- 
istry of H2O. Dashed lines indicate that the reaction is only 
important for temperatures T > 250 K. y denotes the X-ray 
induced FUV photons. The reactions induced by central or in- 
terstellar FUV fields are not included. 



toward the low-mass star-forming region Taurus (TMC 1) of 
Hogerheijde et al. (1997). It is found that = 0.8 x 10 s"' 
fits the observations best. This value is at the lower end of the 
rates found by van der Tak & van Dishoeck ( 2000 1 toward high- 
mass star-forming regions (^„ - 0.61-5.6 x 10 s"'). 
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Fig. 7. Depth dependent temporal evolution of H2O for different X-ray luminosities. The solid line corresponds to the model 
without X-rays. X-ray luminosities are given in ergs s~^ . The T = 100 K mark is at w 10 AU. 



Table 2. Model parameters for a prototypical Class I object. 





0.66 Lq 


Envelope parameters: 




Inner radius ri„ (T = 250 K) 


3.1 AU 


Outer radius (T = 10 K) 


4.3 xlO' AU 


Density at 1000 AU, nCHj) 


8.8 xlO*'cm-3 


Density power-law index, p 


1.6 


Gas mass Me„v {T = 10 K) 


0.034 Mq 


Column density A'h, 


6.9 X 10^- cm-2 



3.2. Results for the Class I model 
3.2.1. X-ray models 

The depth dependent water profiles for f = 5 x lO^yrs to 
t = 10^ yrs and for X-ray luminosities Lx = 10^^-10^** ergs s"' 
are presented in Fig. The inner FUV field is neglected 
in these models. It can be seen that an X-ray luminosity of 
Lx ^ 10^^ ergs s ' destroys the gas-phase H2O within less than 
5000 yrs from initial jcCHjO) ~ 10 ^ down to ^(HzO) ~ 10"^ 
with respect to molecular hydrogen. For an X-ray luminos- 
ity of Lx ~ 10^^ ergs s"' the timescale to destroy water is 
f > 5 X 10^ yrs. In models with Lx < 10^'' ergs s"', water 
has still high abundances (ji:(H20) ~ 10 "*) in the innermost 
part of the envelope. We therefore conclude that X-ray lumi- 
nosities Lx > 10^^ ergs s"' destroy all gas-phase water down 
to fractional abundances of ^(H20) ~ 10"* on timescales 
? 5 X 10"* yrs and Lx > 10^^ ergs s"' on a timescale of a 
few xlO^ yrs. 




Radial Position (AU) 

Fig. 6. Density and thermal structure for a prototypical Class I 
object adopted from the model results of J0rgensen et al. J2002> 
(solid lines). The dashed and dotted lines refer to the gas tem- 
perature for the models including an additional inner source of 
X-ray luminosity (Lx - 10-" ergs s"' and Lx = 10^^ ergs s"', 
respectively). The arrow indicates the distance out to which X- 
rays dominate the H2 ionization rate over cosmic-rays ((a- - 
0.8 X 10-'^ s-i) for Lx = 10^° ergs s"' . 

The water profile for t - lO^'yrs and higher X-ray lu- 
minosities is shown in Fig. |S] The fractional abundances in- 
crease from x(H20) * 10"* to xiH20) ~ 10"^ for increas- 



8 



P. Stauber et al.: Water destruction by X-rays in young stellar objects 



ing X-ray luminosities between Lx - 10^^-10^' ergss and 
;>c(H20) reaches even > 10"^ in the innermost part of the en- 



velope for Lx — 10 ergs s . This effect can be explained by 
the fast electron recombination reaction of H30^ that forms 
H2O and by the increasing gas temperature with increasing 
X-ray flux: The reaction of OH with H2, that drives oxygen 
into H2O, becomes faster than the destruction of water by X- 
rays at higher temperatures, hence more water is produced for 
higher X-ray fluxes. In the region T < 100 K, jc(H20) de- 
creases with increasing X-ray luminosity for f = 10^ yrs, until 
Lx ~ 10^^ ergss Higher X-ray luminosities lead to higher 
water abundances. This is again mainly due to the recombi- 
nation of H30^, which is more abundant in models with high 
X-ray fluxes (see also Sect. Fig. |S] shows that jic(H20) in- 
creases from » 4-5 x 10"^ to a; 2 x 10"^ with increasing X-ray 
luminosity (Lx > 10^^ ergs s"'). At early times, the recombi- 
nation of HjO^ in the models with X-rays leads to even higher 
H2O abundances in the region T < lOOK compared with the 
model without X-rays (Fig.0. 

The OH abundance profiles for diff'erent X-ray luminosi- 
ties at f = 10^ yrs are given in Fig. |9] The OH densities are 
normalized by n(H2) (left) and by «(H20) (right), respectively. 
In the right part of Fig.|9lit can be seen that an X-ray luminos- 
ity Lx = 10^^ ergs s"' leads to n(0H)/«(H20) ^ 10~^ whereas 
X-ray luminosities Lx > 10""^ ergss"' can produce ratios of 
n(0H)/n(H20) * 0.1-1 in the innermost part of the envelope. 
The n(0H)/n(H20) ratio in the hot (L > 100 K) region of the 
envelope in the models without X-rays is only of the order of 
X! 10 X-rays therefore enhance the OH abundances as well 
as the OH to H2O abundance ratio. 

The high densities (hh > 10^ cm""*) and temperatures {T > 
200 K) in the model at r < 10 AU may be appropriate either to 
an inner envelope or to the inner regions of a circumstellar disk. 
However, there is also evidence for cavities in the inner few 
hundred AU (e.g., J0rgensen et al. .2005b) . We have therefore 
also modeled an envelope where Hn is assumed to be at 50 AU 
with a temperature T = lOOK. This corresponds to a source 
with Lboi = 2Lq. The results for t = 10^ yrs are presented in 
Fig.[lO| H2O is destroyed by the X-rays and has abundances 
;i;(H20) < 10"^. At larger radii > 50 AU) the dust will not 
be heated above T > 100 K, i.e., the temperature where water 
evaporates into the gas phase, and the water abundances will 
stay below x(H20) < 10"^ 



3.2.2. FUV models 

The H2O profiles for different inner FUV fields (Go,in = 10^- 
10^) are presented in Fig. ^2 The X-ray emission has been ne- 
glected in these models. As can be seen, inner FUV fields affect 
the water abundances only in the inner 0.2-0.3 AU. This is due 
to the high hydrogen densities (mh ~ 10^ cm ■') in this region 
that absorb the FUV photons very quickly. Although the FUV 
photons from the central source can destroy H2O in the inner- 
most region, even very high FUV fields (Go,in - 10^) cannot 
destroy H2O out to riooK - 10 AU. 

The FUV photons may reach distances in the envelope fur- 
ther away from the source by traveling through the outflow 
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Fig. 9. Fractional abundances of OH (left: x(OH) = 
«(OH)/«(H2), right: x(OH) = «(0H)/n(H20)) for the model 
without X-rays (solid line), for Lx - 10^* ergs s"' (dotted line), 
forLx = 10^" ergss"' (dashed line) and for Lx - 10^^ ergss"' 
(dashed-dotted line). 



cones until they are scattered back into the envelope by dust 
grains (Spaans et al. 119951 j0rgensen|2OO4). However, the in- 
fluence of the scattered photons is restricted to a small region 
around the outflow walls, since the FUV photons are absorbed 
quickly at these high densities and therefore water will not be 
destroyed on large scales. 

Another possibility for the inner FUV field to aff'ect larger 
distances might be a somewhat clumpy envelope. Models of 
clumpy photodissociation regions have shown that the FUV 
photons can travel deeper into the molecular cloud (Meixner 
& Tielens [T993> . At large radii with lower densities, the FUV 
photons primarily heat the gas and not the dust. H2O desorp- 
tion will therefore not be much aff'ected - but H2O gas may be 
photodissociated. Spaans & van Dishoeck <2001l l have shown 
that an inhomogeneous density distribution in molecular clouds 
lowers the column density of H2O by more than an order of 
magnitude compared with the homogeneous cloud model. The 
advent of ALMA or space-based far-infrared interferometers 
might help to investigate the possibility of clumpy envelopes 
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Table 3. Radial column densities of Ht for the Class I source. 
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Fig. 10. Model with ri„ = 50 AU and r(ri„) = lOOK: Depth 
dependent fractional abundances of H2O for different X-ray lu- 
minosities and t - lO'^yrs. The solid line coiTesponds to the 
model without X-rays. X-ray luminosities are given in ergs s"' . 
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Fig. 11. Depth dependent fractional abundances of H2O for dif- 
ferent inner FUV fields and t - 10^ yrs. The solid line corre- 
sponds to the model with Go,in = 0. The X-ray emission is 
neglected. 



in more detail, since a spatial resolution of < 1" is needed to 
resolve the region of interest. 

Models where a central X-ray source and circumstellar 
FUV fields were combined showed similar results for the in- 
fluence of the FUV photons. It is thus concluded that X-rays 
clearly dominate the chemistry of the envelope in comparison 
to FUV fields. 



log(Lx) 




log(Lx) 


A'rad 


log(ergs s-') 




log(ergs s"') 







8.0E11 (2.0E09) 


30 


4.8E14 (4.2E11) 


27 


1.3E12 (2.5E09) 


31 


4.5E15 (3.9E12) 


28 


5.9E12 (6.5E09) 


32 


3.6E16 (2.9E13) 


29 


5.0E13 (4.5E10) 







ionization of CO by He^. Atomic oxygen reaches fractional 
abundances of x(0) ~ 10 in the inner warm (T > 100 K) 
region. The OH^ abundance is dramatically increased by X- 
rays (see also Stauber et al. l2005> . The fractional abundance of 
OH+ is ^ 10 for Lx = lO^^ergss ', x(OH+) ^ 10 for 
Lx = 10^"ergss-i and x(OH+) « 10''°forLx = lO^^e^gss-i 
in the inner region of the envelope. It is mainly produced in the 
ion-molecule reactions of O with , with H2 and OH with 
H+. 

O2 has fractional abundances x(02) ~ 10^^-10 '* in the 
models, much higher than observed. Molecular oxygen is ef- 
ficiently produced in the reaction of atomic oxygen with OH. 
There is an uncertainty in the rate of this reaction for O in the 
ground fine-structure level (^^"2), however, and the rate coef- 
ficient may be lower than that cited by UMIST (Sims et al. 
I2006> . O2 is mainly destroyed in the reaction with H3 and by 
X-ray induced FUV photons. O2 has been searched with SWAS 
(e.g., Melnick 2004i and ODIN (e.g., Pagani et al.|2003l but 
has not been reported in the literature to date. Derived upper 
limits on O2 column densities are usually of the order of » 
a few xlO'^cm"^ or x(02) ~ a few xlO"^. Standard chemi- 
cal models normally overestimate this value by 2-3 orders of 
magnitude (see Melnick 2004 Pagani et al. 2003 for a more 
detailed discussion). Our models are no exception. However, 
the 30" averaged column density in the Class I models are 
A^beam(02) < lO'^ cm"^, which is well within the upper lim- 
its derived by SWAS and ODIN. Oxygen could be frozen out 
on grains in some form at lower temperatures and therefore not 
seen with the big beams of ODIN (9') and SWAS (3.5'x 5.0'). 
Herschel might shed some light on this matter with its relatively 
small beam widths (8-35", depending on frequency). 



Hj is a key molecule in the chemical reaction network 
of water and its abundance is greatly enhanced in the mod- 
els including X-rays. Thus, observations of H3 in absorption at 
a; 3.7 yum (e.g., McCall et al. 1999 1 would allow a direct test of 
the models. The integrated radial column densities A^iad of H, 
are therefore presented for the Class I X-ray models in Table|3l 
The values given in parenthesis are H^ column densities cal- 
culated for the model with ri,, = 50 AU, i.e. the model with a 
protostellar hole. 



3.2.3. Other species 

Depth dependent fractional abundances of the main species in- 
volved in the water reaction network are presented in Fig. 
for t = 10^ yrs. Atomic oxygen is mainly produced by the 
X-ray induced FUV dissociation of O2 and the dissociative 



3.3. Class envelope model 

Table 0] lists the main parameters for a prototypical Class 
envelope model (Schoier et al. I2002> . The crucial differences 
from the modeling point of view between the Class and 
I envelope are the much higher bolometric luminosity, enve- 
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Fig. 12. Depth dependent fractional abundances of relevant species in the chemical H2O reaction network for different X-ray 
luminosities and t - 10^ yrs. The solid line corresponds to the model without X-rays or inner FUV fields. X-ray luminosities are 
given in ergss"'. 



lope mass and column density of the younger object, that lead 
to a different density and temperature distribution. The ini- 
tial chemical abundances are again taken as given in Table [2 
Freeze-out of CO is taken into account for temperatures T < 
20 K. The increased gas temperature due to X-rays has been 
calculated as described in Sect. 13.11 The X-ray luminosity is 
varied between Lx - 10^^-10-'^ ergs s"'. 

The results for H2O are shown in Fig. [O] for t - 10"^ yrs 
and in Fig. [H] for t - 10^ yrs. At t - lO'^yrs, only high X- 
ray luminosities (Lx ^ 10^^ ergs s"') lead to water abundances 
jt:(H20) < 10"^. An X-ray luminosity of Lx > 10^' ergss"' 
destroys water between » 40-80 AU with higher abundances 
WH2O) ~ a few xlQ-^) at « 35 AU and ^ 100 AU. Water 
is not significantly destroyed for X-ray luminosities Lx < 
10^0 ergss Similar conclusions hold for t - 10^ yrs al- 
though the region where ji:(H20) > 10"^ is less than 50 AU 
(T > 250 K). Compared to the Class I models (after appropri- 
ate scaling of the X-ray flux oc Lxr'^), water is less destroyed 
by X-rays in the Class envelope. The reason for this are the 
high densities (ne S 10'' cm"^) in the innermost region of the 
Class source that absorb the X-ray flux. 

4. Discussion 

As discussed in Sect. |2 X-rays reduce the water abundance 
preferably in regions where the gas temperature is L < 250- 
300 K (ji:(H20) « 10"^). At low densities («« < 10^ cm^^), the 
fractional water abundance is less than x 10"^ for X-ray fluxes 
Fx k, 10 ergs s"' cm"- even for high gas temperatures (Fig.|?}. 
For higher densities and temperatures T > 250-300 K, water 



Table 4. Model parameters for a prototypical Class source. 



Lboi 


27 Lq 


Envelope parameters: 




Inner radius n„ 


32 AU 


Outer radius 


8 xlO' AU 


Density at 1000 AU, ;i(H2) 


6.7 xlO'^cm-' 


Density power-law index, p 


1.7 


Mass Mcnv 


5.4 M© 


Column density A'h, 


1.6 X lO^^cm-^ 



has fractional abundances ji(;(H20) > 10 (nearly) independent 
of the X-ray flux. In the following, we discuss these findings in 
the context of protostellar envelopes and other models. 



4.1. Envelope models 

For high-mass YSOs, the observed water abundance can be 
successfully explained within the hot-core scenario in which 
the high temperatures in the inner envelope drive most of the 
oxygen into water and in which most of the water is frozen 
out onto grains in the outer region (e.g.. Doty et al. 120021 
Boonman et al. 2003b;) . Recent interferometer studies of H^O 
by van der Tak et al. J2005> towards the high-mass young stel- 
lar object (YSO) AFGL 2591 confirm the existence of a hot 
(r > lOOK) compact central region with high abundant water 
(x(H20) ^ IQ-*) and a region with lower abundances in the 
colder outer part. Although X-ray models predict the destruc- 
tion of water also in high-mass objects (Stauber et al. i2005l . 



p. Stauber et al.: Water destruction by X-rays in young stellar objects 



11 



10 

o 10" 

CM 

10" 



-A ' 



1 1 


1 1 1 1 1 


1 

Lx =10'° ■ 
= 10 

Lx =10'' ~\ 

Lx =10'' ; 


' \ 

\ 

!.' \ 
■1 \ 
-. \ 


''\ 

1 

\ 


1 \ 
^ \ 


I 
1 


1 t = 10* yrs 

u ~ 


I \ 

\ 


1 




\ 

\, \ 


1 




\ 

\ \._.- 


i 

"^~J..----\ 



/' 


- 

1 1 





50 100 200 

Radial Positioii(AU) 

Fig. 13. Depth dependent fractional abundances of H2O for a 
Class object for different X-ray luminosities and t = 10"* yrs. 
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Fig. 14. Depth dependent fractional abundances of H2O for a 
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The solid line corresponds to the model without X-rays. X-ray 
luminosities are given in ergs s"' . 



these envelopes are likely to have regions with temperatures 



exceeding T = 250 K (van der Tak et al. I2000t . In addition, 
the distance from the central source to the region in the enve- 
lope where the gas temperature is T = lOOK, is much larger in 
high-mass sources compared to low-mass objects due to their 
high bolometric luminosities. The X-ray flux decreases rapidly 
with r^^ and may already be too low to destroy H2O. Water can 
therefore have fractional abundances of the order of 10 in en- 
velopes around high-mass YSOs despite their possible X-ray 
emission. 

Similar jumps in the water abundances are also observed 
toward low-mass Class sources, although the derived abun- 
dances in the inner part are usually much lower compared to 



the massive objects. Ceccarelli et al. (I2000> derived water abun- 
dances of Xin(H20) = 3 X 10"^ for the inner part of the enve- 
lope of the Class source IRAS 16293-2422 and Xout(H20) = 
5 X 10"^ for the outer part from detailed models of ISO obser- 
vations. Maret et al. (2002) found Xin(H20) = 5 x 10"'' for flie 
Class object NGC1333-IRAS4 and XomdiiO) = 5 x 10^1 
Within the spherical envelope models in Sect. 13.31 such jumps 
with relatively low inner water abundances occur for high X- 
ray luminosities Lx - lO-'^-lO"^^ ergs s"' and f - 10^ yrs 
(Fig. I14> . X-rays are therefore a possible explanation for the 
low water abundances found in these objects. 

Class I envelopes may have no or only a small region where 
the temperature exceeds T w 250 K. The inner part of the enve- 
lope may already have been dispersed by outflows or accreted 
onto the star/disk system. Water is therefore easily destroyed by 
X-rays. According to the models in Sect. 13. 2.11 the average wa- 
ter abundance is predicted to be at most x(H20) ~ 10"*' in Class 
I sources with X-ray luminosities Lx > 10^^ ergs s"' . Note that 
the generic Class I model in Sect. l3.ll is considered as an ex- 
treme case where the gas temperature gets T ^ 250 K close 
to the protostar If the inner envelope has cavities of any kind, 
the gas temperature of the envelope will be lower at the inner 
edge of the envelope (^n), since this will be at a larger distance 
and since the dust temperature roughly falls with rT^^'^- Water 
will in that case be destroyed even more readily (Sect. 13.2.11 
Fig. I10> . unless there is some mechanism (e.g., photodesorp- 
tion) that returns water ice to the gas phase. 



4.2. Disks and Outflows 

As discussed in Sect. the inner envelopes of Class and 
I sources are likely to have a much more complex geome- 
try than the spherical envelopes adopted here. If protostel- 
lar cavities or holes (J0rgensen et al. I2005bt are common, 
the emission from H2O and other typical "hot core" species 
may be dominated by the protostellar accretion disk or out- 
flow hot spots rather than by the envelope. Disks are known 
to be present (e.g.. Brown et al. 2QQ0, Looney et al. .2000j 
and can dominate over the hot core in terms of column den- 
sity (J0rgensen et al. '2005c'l. From the chemical point of view, 
the disk and hot core scenario are difficult to distinguish since 
both have high densities and temperatures in their inner region. 
According to the results in Sect. |2j the X-ray flux of a pro- 
tostar with Lx - 10^" ergss"' and an absorbing column den- 
sity of A^H = 10^"* cm"^ at the outer part of a protostellar disk 
(r K 400 AU) is Fx ~ 10"-' ergs s"' cm"^. Water at temper- 
atures below T Si 250 K is therefore destroyed from initial 
jc(H20) ~ 10-^ down to x(H20) * IQ-^-lO ''. High column 
densities and/or temperatures are needed for the water to per- 
sist at x(H20) «i 10 in protostellar disks. For comparison, the 
chemistry disk models including X-rays by Aikawa & Herbst 
( 12001.) have H2O abundances below » 10 at r = 500 AU from 
the central source. In the disk atmosphere models of Glassgold 
et al. (2004 1, water has fractional abundances ji:(H20) < 10"^ 
at 1 AU from the source for vertical column densities of less 
than a few xlO^^ cm"^. 
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Recent results presented by Chandler et al. (I2005> for the 
inner « 1" of IRAS 16293-2422 suggest outflow shocks to 
be the origin for the hot core species. Giannini et al. (2001 1 
and Nisini et al. (2002 1 pointed out the importance of outflows 
for the water abundances in Class and I sources on larger 
scales. The results of Sect. |2 imply that the source of H2O 
emission with ^(HiO) ~ 10""* from outflow regions with gas 
temperatures lOOK < T < 250 K has to be located at a cer- 
tain distance d{Lx) from the X-ray source in order to survive. 
Consider, for example, t - lO'^yrs, typical outflow densities 
(hh - lO'^-lO^ cm"^) and an X-ray absorbing column density 
A'h = 5 X 10^' cm"^. It can then be concluded from Figs.[n- 
13] (Sect. O, that the H2O emitting source would have to be 
at d > 170 AU for Lx = lO^^ergss"', ?A d > 550 AU for 
Lx = lO^'^ergss-i, at £/ > 1700AU for Lx = 10^" ergs s"' 
and at c/ > 5500 AU for Lx = 10^' ergs s '. Larger hydro- 
gen densities and X-ray absorbing column densities, however, 
decrease these distances. The physical quantities of the ob- 
served outflows by Benedettini et al. (.2002J range between 
300K < r < 1400K and 4 x lO^cm^^ < hh < 2 x 10''cm ^ 
It is interesting to see, that the derived water abundances for 
T > 300 K by Benedettini et al. are a factor of k 10 higher 
than those at J = 300 K. This is in good agreement with the 
results in Sect. |2] where the X-ray models show a jump at these 
densities but where the models without X-rays show no jump at 
T - 300 K. X-rays are therefore also predicted to influence the 
water abundances in outflow hot spots. It should be noted that 
FUV fields from the central star may also influence the water 
abundances if the photons can impact the H2O emission region 
relatively unhindered, e.g. through outflow cones. High FUV 
fields tend to destroy water even at high temperatures though 
rSect. l3T2l Stauber et al. l2004l l. 

5. Conclusion 

The gas-phase water abundance is found to be critically de- 
pendent on time, gas temperature, hydrogen density and X-ray 
flux (Sect. Three distinct regimes are identified: In the first 
regime {T < 100 K), the water abundance is ;ic(H20) ~ 10"^- 
10"^ and can be somewhat enhanced or destroyed by X-rays in 
comparison to models without X-rays. In the second regime 
(250 K < T < 100 K), water is released from grains with 
ji:(H20) «i 10""* but quickly reduced by X-rays to fractional 
abundances x(H20) ~ 10 The third regime (T > 250 K) is 
characterized by high water abundances (x(H20) > 10 due 
to the efficient reaction of OH with molecular hydrogen. At low 
densities (nn $10^ cm"^) water is destroyed even at high tem- 
peratures for X-ray fluxes Fx Z 10 ergs s"' cm"^. In general, 
higher gas temperatures and higher hydrogen densities allow 
higher X-ray fluxes for H2O to survive. Water is mainly de- 
stroyed in reactions with the X-ray enhanced species HCO^ 
and H3 . It is also destroyed by internally created FUV photons. 

The envelopes of both Class I (Sect. 13. H and Class 
(Sect. 13.31 1 objects were modeled under the influence of cen- 
tral X-ray and FUV emission. The results for the Class I en- 
velope show that an initial warm gas-phase water abundance 
of x(H20) ~ lO"'* is reduced to jc(H20) ^ 10"^ within ^ 
5000 yrs for X-ray luminosities Lx > 10^^ ergs s ' and within 



=s 5 X 10"* yrs for X-ray luminosities Lx > 10^^ ergs s"' , consis- 
tent with reported upper limits for H2O towards Class I sources. 
Water is also destroyed by X-rays in the Class object, but 
higher X-ray fluxes are needed due to the higher densities in 
the inner region. The influence of a central FUV field is neg- 
ligible in our models, unless the FUV photons can escape to 
larger distances (Sect. llX2l i. 

The current protostellar models are limited to spherical 
symmetry, which is obviously an approximation to the true ge- 
ometry of the inner few hundred AU where cavities and flat- 
tened disks may be present. Nevertheless, the high densities 
and temperatures in the current models are representative of 
such regions. In particular. X-rays are predicted to regulate the 
H2O abundances wherever the gas temperature is T < 250- 
300 K. This will not only be the case for H2O originating from 
envelopes or disks in high-mass star-forming regions, but also 
for H2O produced in outflow hot spots around low-mass ob- 
jects. The presence of cavities or outflow cones will have the 
result that X-rays and/or FUV photons penetrate to larger radii 
in certain directions. 

Future instruments such as Herschel Space Observatory 
will allow detailed studies of the water abundances in both high 
and low-mass YSOs. With its high resolution instrument HIFI 
(de Graauw & Helmich 2001), it will be possible to distinguish 
shock heated H2O from H2O in hot cores or disks from their 
line profiles. Observations of the optically thin isotope H^^O 
will help to sample the envelope, whereas optically thick lines 
may give information about outflow properties. PACS will de- 
lineate the water abundances on 9" scale including maps of 
the outflows. ALMA, on the other hand, will be able to resolve 
the innermost region (< 0.1") of YSO envelopes and/or disks 
in order to study the chemistry and therefore its physical prop- 
erties in molecules other than H2O. Limited maps of Hj^O on 
arcsec scale may be possible though with ALMA under excep- 
tional conditions in the 203 GHz line. High-7 lines of CO will 
reveal the gas temperature in these regions whereas the obser- 
vation of X-ray and FUV tracers (Stauber et al. 2004. 20051 
may help to clarify the high-energy properties of young stellar 
objects. Together, these new facilities and models will reveal 
much about the physical and chemical structure of the inner 
regions of young stellar objects, which are currently poorly un- 
derstood. 
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